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ABSTRACT 

We measure the intrinsic magnitude dispersion of the Galactic bulge red clump 
(RC) using two different methods and arrive at an estimate of o"/ ~ 0.17 mag. We first 
estimate the width of the RC by analyzing in detail a sample of RC measurements toward 
double RC sightlines, which are toward regions of the bulge that are geometrically thin. 
We obtain aj ~ 0.17 mag. We then estimate the value by comparing the luminosity 
function for bulge RC and red giant (RG) stars in V to that in I for a sightline near 
the plane, and we obtain a value of cr/ ~ 0.14 — 0.17 mag. This result has structural 
and evolutionary repercussions. It constrains model estimates of the bar's orientation 
angle derived from RC studies as well as the inferred characteristic length of the bar 
for a given characteristic height. Moreover, the value is too large to account for using 
the predicted evolutionary effects of the known metallicity distribution, implying the 
presence of other population effects. 

Subject headings: Galaxy: bulge 



Introduction 



The red clump (RC) of the Galactic bulge has a vibrant history as an investigative tool to 
enable study of the bulge stellar population, its age, its composition, its distribution, and its dy- 
namics. Its number counts relativ e to the red giant bran ch have been used to either infer a signifi- 
cant intermediate-age popula tion (jPaczvnski et al .11 19941) or mode s tly en hanced helium enrichment 



(|Renzinilll994l : 



Minniti 



19951 ) for the bulge. IStanek et al.l (jl994l . Il997|) used OGLE observations 



of the variation of the dereddened RC mean brightness, brightness dispersion and number counts 
with longitude and latitude to infer that the bulge stars were distributed as bar, with the near side 
oriented toward positive longitudes. The angle between the bar's major axis and the line of sight 
between the Sun a nd the Galactic center, denoted a, was constrained to the range 20° < a < 30°. 
Rattenbury et al.l (j2007al ) obtained similar results using the ~ llx larger OGLE-II survey. 
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Observations of the RC have also been used to trace the r eddening toward the bulge (lUdalski 



2003 
2010 



bum iNishivama et al hood ), to study its dynamics (jRattenburv et al. I l2007bl :] iBabusiaux et al 



De Propris et al.ll201in. and to model the geometry of the bar at large longitudinal separations 



from the G alactic center (jBeniamin et al.ll2005l : iLopez-Corredoira et al.ll2007l : ICabrera-Lavers et al 



20071 . 12OO8I ). It has recently been discovered that the RC is split into two components f or sightlines 
near the bulge minor axis and a t least 5 degrees (~700 pc) removed from the plane (jNataf et al 



2010 



Mc William Sz Zoccalill2010l ). This effect is likely due to the Milky Way's bulge being X-shaped 



at large separations from the plane, a feature predicted by N-body models t o be ubiquitous in b ars 
(|Athanassoulall2005l ) and frequently observed in the bars of other galaxies ([Bureau et al.l 120061 ). 



What allows the RC to be such a powerful diagnostic tool is that it occupies a very nar- 
row range in color-magnitude space with very limited popula tion effects. This is both expected 



from theorv (IGirardi &: Salaris 



2001 



population (Udalski gOOoi r Groeneweg( 
(jPietrzviiski et ahlboiol ). 



:en 



Salaris Sz Girardil l2002l ). observed in the local Hipparcos 



20081 ). and in the stellar populations of other galaxies 



However, while the intrinsic magnitude dispersion of the RC is indeed smaller than that of 
other stellar popul ations, its exact size is not known. In their study of the structure of the bar, 
Stanek et al.l (|l997l ) showed that the axis ratios and orientation angle of the bar are degenerate with 
the assumed value of the intrinsic dispersion. This is due to the fact that the observed brightness 
dispersion is a quadratic sum of the intrinsic dispersion and the geometric dispersion (stars at 
different distances will have different apparent magnitudes), with the geometric dispersion being a 
strong function of the bar structure. A larger assumed value for intrinsic dispersion, and thus a 
lower inferred value for the geometric thickness of the bar, pushes the models to smaller values of 
the axis ratio xq'.zq (a decreased characteristic length for the bar along its major-axis relative to 
its characteristic height) as well as a larger value for a. In this Letter, we estimate the value of the 
intrinsic dispersion of the Galactic bulge RC, that is, the dispersion that would be measured if all 
of the stars were at the same distance and if there was no differential reddening. 

We infer the value of the intrinsic magnitude dispersion using two independent methods and 
estimate it to be aj ~ 0.17 mag. With our first method, we mea sure the brightness dispersion of RC 
stars in the X- wings of the Milky Way bulge using the results of lNataf et al.l (|2010l ) . As the X-wings 
of the bulge are expected from both observations and dynamical models to be substantially thinner 
than the underlying ellipsoid-like structure, the observed dispersion should correspond closely to 
the intrinsic dispersion. With our second method, we compare the observed brightness dispersion 
for the RC in V and / using clump-centric color-magnitude diagrams (CCCMD), a diagnostic tool 
to investigate la rge stellar p o pulatio ns with variable reddening and geometry that is introduced 



and described in iNataf et al.l (2011a 



The structure of this Letter is as follows. We briefly summarize the data used in Section [2l The 
contributing factors to the intrinsic dispersion are estimated from theory in Section [3l A purely 
geometrical estimate using observations of the double clump are summarized in Section [H whereas 



- 3 - 



the method of the CCCMD and the results derived thereof are presented in Section [5j We discuss 
our findings and their implications in Section [6l 



2. Data 



OGLE-III observations were obtained from the 1.3 meter Warsaw Telescope, located at the 
Las Campanas Observatory in Chile. Detailed descri ptions of the in stru mentation, pho t ometr ic 
reductions and astrometric calibrations are available in lUdalskil (j2003l ). and lUdalski et al.l (|2008l ). 



3. The Sources of Intrinsic Brightness Dispersion 



The main contributing factors to the intrinsic dispersion of the bulge RC are expected to 
be the evolution of RC stars during their ~100 Myr lifetime, the effects of dispersion in [Fe/H], 
and that of the dispersion in [a/Fe]. Other factors such as photometric noise, blending due to 
binary contamination, and differential reddening should be no more than a few hundredths of 
a mag each, and will thus not significantly c ontribute in quadrature. Age is expected to be a 
very marginal factor: iGirardi &: SalarisI (|200ll ) estimate that for a ~10 Gyr old , solar-metallicity 
popul ation, dMy/dt and dMi/dt will both be ~0.02 mag Gyr~^. We note that iGirardi k, Salaris 
(|200ll ) predicted an intrinsic dispersion of 0.051-0.107 mag for the bulge RC in I using the metallicity 
constraints available at the time. We refer the reader to their section 5.1. We perform a similar 
procedure here for illustrative purposes, as well as to employ the benefit of more recent data 
constraints. 

The br ightness dispersion a s ingle population would have is estimated using HST observations 
of 47 Tuc (jSarajedini et al.l 120071 ). 47 Tuc is likely to be a single-age, single-metallicity popula- 



tion w ith a small spread (~3%) in initial heli um abundance (Idi Criscienzo et al.ll2010l : iNataf et al 



201 Ibl ). We draw a box around the RC as in iNataf et al.l ( 2011bl ) and find a brightness dispersion 
of 0.065 mag in I and 0.061 mag in V . This value includes the evolution of the RC during its ~100 
Myr lifetime as well as the variation induced by any possible stochastic mass loss during the first 
ascent of the red giant (RG) branch. 



Metallicity dispersion will be a contributing factor to the brightness dispersion. IGirardi &: Salaris 



(|200ll ) predict that for a 10 Gyr population near solar metallicity, dMy/dpe/H] will be 0.60 mag 
dex~^ and (iM//d[Fe/H] will be 0.21 mag dex~^, with Iron-rich stars being fainter. The effect 
of 0.25 dex of enhancement in [a/Fe] is predicted to be 0.15 mag in V and 0.07 mag in I, with 
a-enhanced stars being brighter. The derivatives are thus dMy / d{a /Ye] = —0.6 mag dex~^ and 
dMi/d[a/Ye\ = -0.3 mag dex"!. 



We estimate the metallicity dispersion using the lZoccali et al.l (120081 ) [Fe/H] measurements of 
the upper- RG branch. We find that between the two sightlines h = —4° and h = —6°, the metallicity 
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dispersion changes only slightly, dropping from 0.373 dex to 0.355 dex. Removing stars from the 
distribution with [Fe/H] < —0.80 ~ RG stars that should eventually end up on the RR Lyrae 
instability strip and the horizontal branch but not the RC - reduces the metallicity dispersion 
to 0.300 dex f or b = —4° and 0.278 dex for b = —6°. For [a/Fe], we use the measurements of 
Bensbv et al.l (120101 ) . averaging over [Ca/Fe], [Mg/Fe], and [Si/Fe]. Within the sample, o"[Q,/Fe] — 
0.12 dex, and the value of the correlation coefficient is cor ([a/Fe], [Fe/H]) = —0.86. 

Summing over these factors, we obtain: 



"^/.Predicted 



'^/,single + ^[Fe/K] 



dUi 


2 


L ^2 

r ^[a/Fe] 


dMi 


d[Fe/R]_ 




d[a/Fe\ 



yielding: 



+2cor([a/Fe], [Fe/H])cr[a/Fe]0-[Fe/H] 



a"7Predicted = 0.117mag 



dMi (Mi 
d[a/Fe] d[Fe/R] 



(1) 



A similar calculation for fXy^Predicted yields a value of 0.256 mag, for a predicted intrinsic brightness 
dispersion ratio of av/cri ~ 2.18. 

Our theoretical estimate is done under the assumption that the bulk of Galactic bulge stars 
formed in a very small time-frame without enhanced helium-eri r ichment. The stars t hat are rich 
in [Fe/H] also have lower values of [a/Fe] ([Bensbv et al.l l20ld : iGonzalez et al.ll201ll ). as such a 



strong prediction of this single-age, simple-helium hypothesis is t hat the most rn etal-rich RC stars 
should on average be the dimmest. This trend is observed locally. lUdalskil (|2000l ). in his analysis of 
Hipparcos RC stars with metallicity measurements, found that RC stars with [Fe/H]~0.0 were 0.08 
mag fainter than RC stars with [Fe/H]~ —0.36. If the most metal-rich stars are either substantially 
younger or more helium-enriched, then that correlation may be washed away or even reversed. 



We note the main reason why we obtain a larger estimate than lGirardi &: Salarid (|200ll ). It is 
primarily due to the the last term in Equation (1), the anti-correlation between [a/Fe] and [Fe/H]. 
The stars with the highest [Fe/H] also have the lowest [a/Fe], making them even fainter than 
expected relative to the stars with lower [Fe/H]. 



4. Intrinsic Brightness Dispersion from Observations of the Double Clump 



The double RC of the Galactic bulge (iNataf et al.ll2010l : lMcWilham fc Zoccalill20ld ). in which 
two RCs with equal {V — I) and (J — K) color but with a 0.5 mag separation in brightness coexist, 
is a phenomeno n observed toward si ghtlines (|^| ^ 3.0, \b\ ^ 5.0). X-shape d bulges are ubiquit ous in 
N-body models (|Athanassoulall2005l ) and in observations of field galaxies ([Bureau et al.ll2006lV The 



two RCs show a significant mean difference in longitudinal proper motion (|Mc William &: Zoccali 



20 id ) and are accompanied by two RGBBs of similar brightness and number counts ([Nataf et al 
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201 lal ). In light of these factors, we are confident in making the assumption that this is a geometric 
effect. 

Studies of X-shaped bulges in other galaxies indicate that the geometric al dispersion toward 
these sightlines should be very small. The examples shown in Figure 1 of (jBureau et al.l l2006l ) 
have X-wings with total length typically ~l/4 that of the bar's half-length. For a ~2.5 kpc Milky 
War ( Rattenburv et al.ll2007al ). that implies a total length of ~600 pc. Dividing by ^/l2 gives a 
maximum dispersion of 175 pc, or ~0.05 mag for a structure at 8 kpc. We note that this is an 
upper-bound as the value of \/T2 is only correct if the stars are uniformly distributed within the 
X-wings, an unlikely propostion. If they are more heavily concentrated toward the center, the 
ratio of geometrical dispersion to total physical extent will be smaller, and the geometric dispersion 
consequentially reduced. 



We use a sample of 591 double RC measurements used as part of the study of iNataf et al 



(j2010l ). We fit for the KG + double RC stellar population using the equation: 



N{m) = ^exp 



B{m-mRc,i) 



+ 



Nrc: 



2vr cr/jc, 



exp 



[m 



mRc,i) 



2n 



RCA 



+ 



Nrc- 



exp 



(m - ruRc^f 
RC,2 



(2) 

where the parameters are as described in lNataf et al.l ( 201Cl ). For two Gaussians placed near each 
other, a fit will often have a reduced jf tj^e normalization and dispersion of one Gaussian is 
increased at the expense of that of the other Gaussian. We correct for this by imposing the 
geometrically-motivated prior that the two Gaussians should have the same brightness dispersion. 
We show a representative CMD in Figure [H 

We obtain a mean dispersion of 0.174 mag in / for the two RCs. We plot the distribution of 
magnitude dispersions toward these sightlines in Figure [21 We note that the standard deviation 
in the mean value of the dispersion obtained for these fields, 0.046 mag, is roughly equal to the 
average of the standard deviations for the width of the RCs in each of the Markov Chains used to 
obtain the best fit, 0.041 mag. This implies that all or nearly all of the variance seen in Figured] 
is due to Poisson noise rather than structural changes. 

Subtracting in quadrature the expected geometric dispersion of ~0.05 mag from the obtained 
mean of 0.174 mag yields 0.167 mag. Therefore, our analysis of the double RC sightlines yields an 
intrinsic RC brightness dispersion o"/ ~ 0.17 mag. 



1 1.5 2 2.5 

V-I 



Fig. 1. — The double RC toward the sighthne (l,b) = (0.0, —6.3). The two peaks in the magnitude 
distribution of RG+RC stars are clearly discernible at / = 14.79 mag and / = 15.27 mag. The 
right panel shows the magnitude distribution of the RG+RC stars in / for the subset of the stars 
in the left panel that are within the black shaded contour lines selecting the RG branch. 
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Fig. 2. — Distribution of Gaussian widths urc obtained for a sample of 591 sightlines toward the 
Milky Way's X-wings. 



5. Measuring the Intrinsic Brightness Dispersion Using Clump-Centric Color 

Magnitude Diagrams in V and / 

We use CCCMDs to study the RC for OGLE-III point sources that are in OGLE-III subfields 
centered within 0.5 degrees of {I, b) = (0, —2). The construction of the CCCMDs, in which the color 
and magnitu de of each s t ar is ta ken relative to that of the nearest measured RC centroid, is fully 
described in iNataf et alj (|2011al ). Of the available sightlines in the OGLE-III bulge photometric 
catalog, {l,b) = (0,-2) is among the ones that is closest to the plane, and will as such have the 
least amount of geometrically induced brightness dispersion for single-clump fields. The combined 
brightness+intrinsic dispersion is ~0.22 mag, compared to ~0.25 mag near Baade's window. This 
narrower dispersion allows a cleaner separation between the RC and the re d giant branch bump 
(RGBB), reducing the degeneracy between the two. As in lNataf et al. ( 2011a ). we fit for the various 
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components of the stellar population as follows: 



N{m) = Aexp 



B{I - Irc) 



, Nrc 
H — 7^ exp 



2'KaRc 



{I 



Irc? 



'^^RC 



Nrgbb 
H — 7^ exp 



'^'^(^RGBB 



{I - IrgbbY 

^^RGBB 



(3) 

where vl, B are the parameters to the underlying RG branch, Nrc, Irc, and grc are the normal- 
ization, mean brightness and brightness dispersion for the Gaussian used to model the RC; and 
Nrcbb, Irgbb, and ctrcbb are the analogous parameters for the RGBB. About 230,000 sources 
are used to fit the combined RC+RG+RGBB population. We obtain aRc = (0.224 ± 0.002) mag 
in I, with the quoted error being a statistical error. The resulting fit is shown in the top panel of 
Figure [3l 



CT„„ = (0.2^4 ±0.002) mag 
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cTrc = (0.3i36±'o.o63) mag 




V - V 
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Fig. 3.— TOP: Fit for the magnitude distribution in / of ~230,000 OGLE-III RG branch point 
sources near (/, b) = (0, —2). BOTTOM: Same as top, but the magnitudes used are in V. 
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We repeat the process in V though we constrain the RGBB parameters as they are less stable 
when the RC is broadened. The population ratio fsump = ^n%c^ fixed to have the same value 
as obtained when fitting in /, the magnitude difference is fixed to be the same (an assumption of 
identical mean {V — I) color for the RGBB and RC), and (Trgbb is fixed to 0.3 mag. The resulting 
fit is shown in the bottom panel of Figure [3l We obtain cry = (0.336 ± 0.003) mag. We note that 
the value obtained here is very weakly sensitive to the reasonable priors we have imposed on the 
RGBB in V . Removing the constraint that the RGBB and RC to have the same mean color raises 
ay from 0.336 to 0.343 mag, whereas floating (Trgbb as a free parameter in V raises cry from 0.336 
to 0.347 mag. 

A few results are immediately apparent from Figure [3j First, the RC is much more well-defined 
in I. There is no genuine prospect, for example, of ever measuring the parameters of the Galactic 
bulge RGBB in V . The secondary over density in simply not discernible in V , and its parameters 
must be imposed. 

We comment on the validity of the assumption that the brightness of the RC is distributed as 
a Gaussian. It is known from globular cluster studies that the RC of a single-age, single-metallicity, 
single-distance population is intrinsically non-Gaussian - there's a steep rise in number counts at 
the location of the zero-age horizontal branch, followed by a slow gradual decline. However, these 
sightlines have additional noise due to the broad metallicity variation a nd geometrical th i cknes s 



of the Galactic bulge. Indeed, the simulated bulge RCs in Figure 14 of iGirardi &: Salarid (120011 ) 
have very Gaussian-looking brightness distributions, and that is without the additional effect of 
geometrical broadening. 

The intrinsic magnitude dispersion cr/ can be estimated from the combined geometric and 
intrinsic magnitude dispersion in / and V . We assume that the RC in / and V will have the same 
contribution from geometric dispersion, but different contributions due to population effects: 

4c,/ = ^G + '^^i (4) 

2 2,2 fr-\ 

Where ac is the geometric dispersion, aj is the intrinsic dispersion in /, and analogously for cry. 
Substituting in aj = 0.224 and cry = 0.336 and taking the difference of the equations yields: 

al-a] = 0.0672 (6) 

The value of aj is then a monotonic function of the ratio cry / aj . 

Adopting the expected value of oyjoi = 2.18 derived in Section [3] yields cr/ = 0.14 mag. A 
consideration of the uncertainties in this estimate - and thus of the uncertainty in uy/cr/ - is in 
order. Of the effects listed, the only one that in and of itself produces a ratio less than oy joj = 2.00 
(for which o"/ = 0.15) is that of the luminosity evolution of a RC star during its lifetime. It is ~0.06 
mag for 47 Tuc. If the value were doubled, to 0.12 mag, the ratio would be lowered to ~1.8, yielding 
cr/ = 0.17. Alternatively, using the value of cr/ = 0.17 measured in Section |4] yields ay /aj = 1.8. 
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We therefore adopt aj ~ 0.14-0.17 mag as the estimate obtained by comparing the dispersions 
in V and /. 



5.1. On the Assumption of Minimal Differential Reddening 



At the request of the anonymous referee, we discuss the assumption of minimal differential 
reddening within the sightlines used to construct the CCCMD used in section 5. 

We first note that the sightlines used to measure the RC centroids before building the CCCMD 
are no more than 3' apart, which means that any point source kept will be within ~2' of a measured 
RC centroid . This maxi i nizes a ny differential reddening to that leftover on ~2' scales. Second, as 
discussed in lNataf et al.l (|2011al ). when constructing our clump-centric photometry files we directly 
visually inspected each sightline and discarded those with very high differential reddening, for 
example those where the RC would have a slanted profile on the CMD. 



dA,/dE(V-I) = 0.99 




Fig. 4. — The color and magnitude of 304 RC centroids measured within 0.5 degrees of = 
(0,-2). 

The 304 RC centroids measured within 0.5 degrees of the sightline {l,b) = (0,-2) have a 
standard deviation in coordinate space of ~15', and a standard deviation in reddening of 0.4 mag 
in E{V — I). We show the scatter plot in Figured! This implies that even if the most egregious 
cases of differential reddening were not directly removed, the expected differential reddening per 
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sightline would be ~0.05 mag in E(V-I) and thus the same in / due to the reddening law toward 
these sightlines, Ai/E{V — I) = 0.99. As the intrinsic dispersion estimated in Section 5 is o"/ ~ 
0.14-0.17 mag, differential reddening will not significantly contribute in quadrature. 



6. Conclusion 



In this Letter, we obtain similar values for the intrinsic brightness dispersion of the Galactic 
bulge RC using two different methods, by analyzing the fits to double RCs toward (—3.5 ^ I ^ 
1.5, \h\ ^ 5.0), and by comparing CCCMDs toward (/, h) ?a (0, —2) . We respectively obtain values 
of a I 0.17 and ai 0.14 — 0.17 mag using the two di fferent methods. Our result has consequences 
for studies of Galactic structure. IStanek et al.l (jl997l ) found that decreased values of the intrinsic 
magnitude dispersion of the RC would yield decreased values of the bar angle a and the bar 
structural parameter xq:zq. 



m 



These value C7f ~ 0.17 mag is larger than the theoretical estimates of ~0. 051-0. 107 mag found 
Girardi &: SalarisI ()200ll ). It is also higher than our theoretical estimate of ~0.117 mag. The 



difference may be due to errors in the theory, binary evolution effects, or correlated variations 
between the input parameters we considered with age and helium-enrichment. 
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